Distillation of Entanglement between Distant Systems by Repeated Measurements on 

Entanglement Mediator 



G. Compagno/'Q A. Messina, 1 [j] Hiromichi Nakazato, 2 '[| A. Napoli, 1 '!! Makoto Unoki, 2 and Kazuya Yuasa 2, 

1 INFM, MURST and Dipartimento di Scienze Fisiche ed Astronomiche 
dell'Universita di Palermo, Via Archirafi 36, 90123 Palermo, Italy 
2 Department of Physics, Waseda University, Tokyo 169-8555, Japan 

(July 31, 2004) 

A recently proposed purification method, in which the Zeno-like measurements of a subsystem 
can bring about a distillation of another subsystem in interaction with the former, is utilized to yield 
entangled states between distant systems. It is shown that the measurements of a two-level system 
locally interacting with other two spatially separated not coupled subsystems, can distill entangled 
states from the latter irrespectively of the initial states of the two subsystems. 
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I. INTRODUCTION 

One of the key technologies for quantum information 
and computation is purification/distillation of quantum 
states pi 12- Particular pure states, such as entangled 
states, often play significant roles there, but it is not 
easy to find such "clean" states in nature. It is therefore 
required to prepare them out of mixed states; otherwise, 
we can not carry out any interesting ideas of quantum 
information and computation. 

A new purification mechanism has recently been pro- 
posed 0. It is shown that repeated measurements on a 
system, say A, result in a purification of another system, 
say B, in interaction with A |4|. That is, the state of B 
is driven to a pure state irrespectively of its (generally 
mixed) initial state, if certain conditions are satisfied. 
Remarkably, if appropriate adjustment of the relevant 
parameters is possible, the maximal yield, which is pre- 
scribed by the initial mixed state of B and its target pure 
state, can be attained, while keeping the maximal fidelity, 
by a finite number of measurements on A (an "optimal 
purification" ) . This constitutes a remarkable contrast to 
the standard purification protocol P, Q, in which it is 
generally difficult to realize both a non-vanishing yield 
and the maximal fidelity at the same time. 

Since an entangled state is one of the pure states of 
two quantum systems, say A and B, one can think of 
a possibility of extracting the entangled state between 
A and B by repeatedly performing measurements on X 
which interacts with both A and B. This possibility has 
already been pointed out and explored to show that 
one of the Bell states can be extracted when this mecha- 
nism is applied to a three-qubit system, where qubits A 
and B always interact with the other qubit X on which 
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one and the same measurement is repeatedly performed. 
Notice that in this case the two systems A and B are not 
spatially separated, because they are supposed to locally 
interact with X. On the other hand, it is often required, 
e.g., in the ideas of quantum teleportation and commu- 
nication 0j j to establish an entanglement between two 
quantum systems that are located at or at least can be 
sent, without loosing the entanglement, to distant places. 
In this respect, it would be worth remembering that in- 
teresting ideas of generating an entanglement between 
two cavities |8( and of transferring an entanglement be- 
tween two modes in a cavity to that between other two 
modes in different cavities Q have been proposed. In 
the former a two-level atom is sent to interact succes- 
sively with the two cavities resulting in the generation of 
an entanglement between the two, and in the latter the 
entanglement is shown to be transferred by a two-level 
atom which passes through the two cavities and interacts 
with the relevant cavity modes. In this paper, this kind 
of successive interaction with two quantum systems is in- 
corporated within the framework of the new purification 
mechanism 0, 0, |^ to show that an entanglement can be 
established between the states of the two systems spa- 
tially separated (or that can be separated). Notice that 
the entangled state is distilled from an arbitrary initial 
state that is in general mixed, while in the generation of 
entanglement in Ref. the initial state should be pre- 
pared in an appropriate pure state and in the transfer 
of entanglement in Ref. |9( the state is assumed to be 
initially entangled. 

After a brief review of the new purification mechanism 
in Sec.m a scheme of successive interaction is introduced 
in a three-qubit system, A+B+X, in which system X is 
assumed to interact first with system A and then with 
B, in Sec. II I II System X is prepared in an initial pure 
state and is measured after it has interacted with A and 
B. Then only those events in which system X is found in 
the initial state are kept. This process will be repeated 
many times. It is shown that an optimal entanglement 
purification is actually realizable for a particular choice 
of interaction and by properly adjusting interaction times 
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and strengths between A and X and B and X. In Sec. II VI 
another example of an entanglement purification is ex- 
amined in a physical system where a two-level atom X is 
injected successively to the two cavities A and B "back 
and forth," interacts with their cavity modes under the 
rotating-wave approximation and the state of X is re- 
peatedly measured in a prescribed way. It is explicitly 
shown that, under certain conditions, a particular entan- 
gled state between the lowest two modes of each cavity is 
extracted, irrespectively of the initial cavity states. Fi- 
nally we summarize the results obtained and give future 
perspectives in Sec. [S 



II. PURIFICATION VIA REPEATED 
MEASUREMENTS 

Let the total system consist of two parts, system A 
and system B, and the dynamics be described by the 
total Hamiltonian 



H = H A + H B + H h 



(2.1) 



where H mt stands for the interaction between the two 
(sub)systems. We initially prepare the system in a prod- 
uct state 



po = |0>(0| ® p B (0) 



(2.2) 



at t = 0. Notice that system B can be in an arbitrary 
mixed state pb(0). We perform measurements on A at 
regular intervals r to confirm that it is still in the state | <t>) 
[4| , while the total system A+B during the time r evolves 
unitarily in terms of the total Hamiltonian H. Since the 
measurement is performed only on system A, the action 
of such a (projective, for simplicity) measurement can be 
conveniently described by the following projection oper- 
ator 



(2.3) 



Thus the state of system A is set back to \4>) every af- 
ter r, while that of B just evolves dynamically on the 
basis of the total Hamiltonian H . We repeat the same 
measurement, represented by (|2.3|) . N times and collect 
only those events in which system A has been found in 
the state \<p) consecutively N times; other events are dis- 
carded. The state of system B is then described by the 
density matrix 



P ( b\n) = (V,(r)) N PB (Q)(vJ(r)) N /P^(N) 
where 

is an operator acting on B and 

P^(N) = Tr\(Oe- lHT 0) N P0 {Oe lHT O) N 



(2.4) 



(2.5) 



= Tr, 



^(t))%b(0)(V/(t))"1 (2.6) 



is the probability for these events to occur (yield). This 
normalization factor appearing in 12.4fl reflects the fact 
that only the right outcomes are collected in this process. 

In order to examine the asymptotic state of system B 
for large N, consider, assuming its existence, the spectral 
decomposition of the operator V$ (r) , which is not hermi- 

tian, Vtfjlr) ^ VUt). We therefore need to set up both 
the right- and left-eigenvalue problems 

V^(r)|u„) = \ n \un), (« n |V^(r) = A„(v„|. (2.7) 

The eigenvalue A„ is complex valued in general, but its 
absolute value is bounded |5j 



< IAJ < 1. 



(2.8) 



This reflects the unitarity of the time evolution operator 
e~ T . These eigenvectors are assumed to form a com- 
plete orthonormal set in the following sense 

22 1 u n )(v n | = 1 B , (v n \u m ) = S nm . (2.9) 



(We normalize \u n ) as (u n \u n ) — 1, while the norm of 
(v n | has been fixed by the above relations and is not nec- 
essarily unity.) The operator (r) itself is now expanded 
in terms of these eigenvectors 



(2.10) 



It is now easy to see that the Nth power of this operator 
is expressed as 

(V4r)) N = y2^\u n )(v n \ (2.11) 



and therefore it is dominated by a single term for large 
N 



N largo N N . . . . 

> A Fo)(«o|, 



(2.12) 



provided the largest (in magnitude) eigenvalue Ao is dis- 
crete, nondegenerate and unique. If these conditions are 
satisfied, the density operator of system B is driven to a 
pure state 



( T )/AT\ larg ° N I W I 

Pb ( n ) > K)(u | 



with the probability 

pM(tf) largc jV , 



X r (v o \pb(0)\v ) 



(2.13) 



(2.14) 



The pure state \uo), which is nothing but the right- 
eigenvector of the operator V$ (r) belonging to the largest 
(in magnitude) eigenvalue Ao, is thus distilled in system 
B. This is the purification scheme proposed in 

A few comments are in order. First, the final pure 
state |uo) toward which system B is to be driven is de- 
pendent on the choice of the state \<p) on which system 
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A is projected every after measurement, on the measure- 
ment interval r and the Hamiltonian H, but does not 
depend on the initial state of system B at all. In this 
sense, the purification is accomplished irrespectively of 
the initial (mixed) state /?b(0). Second, as is clear in 
the above exposition, what is crucial in this purification 
scheme is the repetition of one and the same measure- 
ment (more appropriately, spectral decomposition) and 
the measurement interval r need not be very small 0. 
It instead remains an adjustable parameter. Third, if 
we can make other eigenvalues than Ao much smaller in 
magnitude 



|A n /A |«l for n ^ 0, 



(2.15) 



by adjusting parameters, we will need fewer steps (i.e., 
smaller N) to purify system B. 

It is now evident that the purification can be made 
optimal, if the conditions (|2.15(l and 



|Ao| = 1 



(2.16) 



are satisfied. This condition l|2.16[) assures that we can 
repeat as many measurements as we wish without run- 
ning the risk of losing the yield (N) in order to make 
the fidelity to the target state \uo), 



F&(N) =Ttb \ P <£\n)\u )(u \\ , 



(2.17) 



higher. Actually, the yield p( T \N) decays like 

pM(JV) =^2\%\* m N (v n \p B (0)\v m )(u m \u n ) 



|A | 2A >o|pb(0)N> 



(2.18) 



and the condition (|2.1tj[l can bring us with the non- 
vanishing yield (wo|pb(0)|«o) even in the N — > oo limit. 
Therefore the condition (|2.16() makes the two (sometimes 
not compatible) demands, i.e., higher fidelity and non- 
vanishing yield, achievable, with fewer steps when the 
condition (|2.15|l is met. In this sense, the purification is 
considered to be optimal. 

It would be desirable if an optimal purification can be 
realized by an appropriate choice of the state \<j>) and/or 
tuning of the measurement interval r and parameters in 
a given Hamiltonian. A few simple systems have already 
been examined 0, 0, |(| to show how such optimal purifi- 
cations are made possible. 



III. ENTANGLEMENT DISTILLATION IN A 
TWO-QUBIT SYSTEM A+B BY AN 
ENTANGLEMENT MEDIATOR X 

Since an entangled state of system A+B is one of the 
pure states, there is a possibility that we apply the pu- 
rification mechanism described above to distill the initial 
(generally mixed) state of A+B to a desired entangled 



state. This possibility has already been pointed out and 
it has been explicitly demonstrated that one of the Bell 
states | V I'~) of the two qubit systems A+B can actually 
be extracted if we repeatedly measure one and the same 
state of another qubit system X, the interactions of which 
are symmetrical with respect to A and B, resulting with 
the maximal yield 0, HI • One of the limitations of this 
model is that the entanglement can be established only 
when the two systems A and B locally interact with the 
same system X at the same time and therefore it does 
not seem to allow to establish an entanglement between 
two systems that are spatially separated. A new scheme 
is certainly needed. 

When the two systems A and B are spatially separated, 
their local interactions with the other system X can not 
take place simultaneously. This means that the inter- 
actions are considered to become effective one by one, 
i.e., system X first interacts with, say system A and then 
with B 0, [|| . This kind of process can be conveniently 
described by a time-dependent total Hamiltonian. In this 
section, a total system composed of three qubits (or three 
spin 1/2 systems), A, B, and X, is considered. The two 
qubits A and B interact with the other qubit X succes- 
sively, the state of X, initially prepared in a particular 
state, say in up state, is measured after the interactions 
with A and B, and only those events in which the state of 
X is found in up state will be retained and other events 
are just discarded. The process will be repeated many 
times and we are interested in the resulting state of A+B. 

Assume that the three-qubit system A+B+X is de- 
scribed by a time-dependent total Hamiltonian H(t). 
Qubit X, which is initially prepared in up state |Tx)> i s 
first brought to interaction with qubit A for time interval 
<A- The Hamiltonian in this period is 



H{t) = H Q 



H XA- 



(3.1a) 



Next, after a free time evolution under the free Hamil- 
tonian Hq for time duration ta, qubit X interacts with 
another qubit B, which has no direct interaction with A, 
for ie- The Hamiltonian for this period reads 



H(t) = H Q 



H XB- 



(3.1b) 



After another free time evolution for tb , the state of X is 
measured and only those cases in which qubit X is found 
in its initial up state || x ) shall be retained. The whole 
process, i.e., 

interaction between X and A for t\ 
— > free evolution for ta 
— > interaction between X and B for £b 
— ► free evolution for tb 

— > projection to |T X ) (3-2) 

shall be repeated many times and we are interested in 
the final state of A+B. See Fig. 1. 

In order to describe the above process explicitly and to 
show the possibility of entanglement distillation by this 
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A 



0-0^ 



ITx) 



FIG. 1: Qubit X, prepared in |T X ), is brought to interaction 
with qubits A and B successively and then its state is mea- 
sured. If it is found in |tx)> the whole process is repeated 
again; other events are discarded. 



process, we consider the following free and interaction 
Hamiltonians 



(1+4 A) ) + ^(1 + 4 B) ) + ^(1 + 4 X) ), 



H - 

rjl (X) (A) 



H XB 



(X) (B) 

9b<?\ o\ , 



(3.3a) 
(3.3b) 



(A) 

where a\ s are Pauli matrices acting on the Hilbcrt 
space of system A and so on, and g A and are real 
(assumed, for definiteness, to be positive) coupling con- 
stants. The free Hamiltonians for A, B, and X are as- 
sumed to be the same for simplicity and are characterized 
solely by the common energy gap u>. 

It is clear that the relevant evolution operator for the 
whole process (|3.2[1 is given by 



iH T-B -i(H +H' XB )t-B 



x e - iHoTA e- i< - Ho+H ^ )tA ^ x ). (3.4) 



It is an elementary task to evaluate this operator, since 
each factor on the right-hand side is easily evaluated in 
terms of the eigenstates of the Hamiltonian in the expo- 
nent. Indeed, we have 



i(ffo+-ffxA)*A 



It, 



_ -iwt A (l+CT^ ')/2-itvtA 



and similarly 

/j | e - i ( H o+-Hi B ) t B 



-iujt B ( l+<? { 3 ' ) /2-iwt B 



I Tx> { (cos <PA ~ i sin ip A cos 20 A ) | Ta> <T a I + cos(g A t a) 1 1 a) (I a I } 

+ Ux){-* sin ^A Sm26> A || A )(T A | ~ * sin (s A £ A )|T A )(l A |} 

(t x |{(cos(p B - isin<^B cos20 b )|Tb)(TbI + cos(5 B i B )|l B ) (IbI> 
+ (lx\{- isin( PB sin26» B |t B )(lBl - *sui( 5 b£b)Ub)(TbI} 
I 



(3.5a) 



(3.5b) 



where the angles <^a(b) and $a(b) are denned as 



<PA(B) = t A (B) \Ju 2 + 9 A (By 
9A(b) 



sin26» A ( B ) = 

C0S26» A (B) = 



^ 2 + 9 2 MB ) 



9 A (B) 



(3.6a) 
(3.6b) 

(3.6c) 



Let us introduce a parity operator V = <7g 03 whose 
eigenvalues +1 and —1 single out two subspaces of the 
product Hilbert space TL A ®H B invariant under the action 
of the operator V. The two states |TaTb) and | Ia-Lb) 
generate the even parity subspace and the following 2x2 
matrix M. with its elements 

Mu = e" lw(tA+2rA+tB+2rB) (cos^ A -isin(^ A cos26l A ) 
x (cos^b — isin^B cos2#b), (3.7a) 



M 12 
M21 



sin (p a sin 26 'a sin (#b t b ) , 



(3.7b) 



-M*b+2t b ) s i n ( 3A t A ) s i n(fiB s - m 29 B , (3.7c) 



M22 = cos(g A t A ) cos(g B t B ) 



(3.7d) 



allows to completely characterize the action of V in this 
subspace as follows 



V 



ITaTb) 
IIaIb). 



g — iu(tA + TA + tB+TB)^/f 



ITaTb) 
IIaIb). 



(3.8) 



We proceed in the same way for the odd parity subspace 
spanned by the states |TaIb) and |IaTb)- To this end we 
define the 2x2 matrix M with its elements 

Mi = e~^ (2rA+tB) (cos ip A - ismip A co82e A ) cos(g B t B ), 

(3.9a) 

A/12 = — sin if A sin 20 A sin ip B sin 20b , (3 . 9b) 

AA 21 = _ e — (*A+2r A+fB ) sin ( 5AU ) si n( 5B te), (3.9c) 

M22 = e -^(* A+2TA ) cos( ffA i A )(cos^ B - isin</>BCOs20 B ), 

(3.9d) 

so that the action of V is represented as 



V 



ITaIb) 

UaTb) 



j(tA+t B +2T B )j\j- 



IT a 4b> 
IUTb) 



(3.10) 
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In order to show explicitly that the process l|3.2[l with 
the particular choice of interaction (|3.3b|l admits an en- 
tanglement distillation for qubit-system A+B, it turns 
out to be enough to consider a much simpler case. Let 
us treat systems A and B symmetrically, except for the 



ordering of their interactions with system X. We choose 
the same parameters for A and B, i.e., <?a = 9b = 9, 
t A = t B = t and r A = r B = r (lpa(B) — > <p and 
^A(B) ~~ * F° r the parity-odd states, the matrix M 
now is simplified to be 



[e 4 "(*+ 2t ) cos(gt)(cos tp — i sin tp cos 20) 
_ e -M"(frH-)ain 2 (fl*) 



— sin 2 y> sin 2 20 ^ 

e -tu{t+2r) cos (gt) ( cos _ j sm ^5 cos 26) J 



(3.11) 



It is easy to find the condition under which an entangled 
state of the form 

|*)^-i=(|T A i B }+ e * x || A TB», (3-12) 

where x is a real parameter, is an eigenstate of this matrix 
N (and therefore, of the operator V). A straightforward 
calculation shows that if the parameters g, t and r are 
so chosen that the following relation 



cos tp — i sin tp cos 26 = —e luJT cos(gt), 



(3.13) 



is satisfied, the state \ty) with \ = u>(t + r) is in- 
deed an eigenstate of V belonging to the eigenvalue 
Ao = — e ~ 3M (*+ r ). [There is another possibility of op- 
timal distillation of the above entangled state (|3.12|) , but 
with a different \, i.e., \ = UJ ('t + t) + n. This case is 
realized under the condition (|3.13l) with the replacement 
lot — > lot + 7r; the corresponding eigenvalue is also given 
by the shifted one, i.e., e - 3 ^(*+ T ).] 

Notice that we are not allowed to set cos(gt) s\n(gt) = 
because it would result in a degenerate (in magnitude) 
eigenvalue. Observe that we have essentially two condi- 
tions l|3.13|l . while we have three independent combina- 
tion of parameters gt 1 Lot and lot. We, therefore, have 
a possibility of an optimal distillation of the entangled 
state \^>) , if the magnitudes of the other eigenvalues of V 
are made smaller than unity. The remaining eigenvalue 
of M under the conditions in i|3.13|l reads 



JT cos(gt) (cos tp— i sin tp cos 26) +sin 2 (gt) 
e- iw ^ [- cos 2 (gt) + sin 2 (gt)], (3.14) 



and its absolute value can not be made unity when 
cos(gt) sin(gt) ^ 0. On the other hand, matrix M. is 
expressed as 



-2iu{t+T) 



M = 



cos 2 (gt) =Fe" l "*sin 2 (gi) N 



^- Fe - lW (t+2r) gin 2 ^ CQS 2 (gQ 



(3.15) 



and the absolute values of the eigenvalues of this ma- 
trix can not reach unity if cos co(t + r) ±1 and 
cos(gt) s'm(gt) ^ 0. This means that, under these con- 
ditions on gt, tot and lot satisfying the relation l|3.13|l . an 



optimal purification (i.e., distillation) of the entangled 
state (|TaIb) + e M( * +T) || A TB>)/V2 is possible. It is in 
fact easily shown that the left-eigenstate of V belonging 
to the eigenvalue Aq = — e ~ 3 ' iW ( t + T ) is expressed as 



M>l 



v/2 



((TaIbI + ^UaTbI), (*|*>= 1, (3-16) 



and therefore the yield P^(N) approaches asymptoti- 
cally, as N becomes large, a finite value 

P W (oo) = i[<T A | B |+e-^aAt B | 

xpAB(0)[|T A 4B>+e iX |lATB>]. (3-17) 

This is nothing but the probability of finding the target 
entangled state |*) = (|Ta!b) + e ix |J. A T B ))/V2 in the 
initial state pab(0). 



IV. ENTANGLEMENT DISTILLATION OF 
CAVITY MODES 

In the previous section, the possibility of realizing an 
entanglement distillation is demonstrated for the three- 
qubit system A+B+X. The particular form of interac- 
tion l|3.3b(l is shown to be suitable for this purpose follow- 
ing the procedure (|3.2|) . In this section, another applica- 
tion of the purification mechanism 0, [3, |fj is explored in 
a system composed of a two-level system (e.g., an atom) 
interacting with two single-mode cavities. The two cav- 
ities may be located at spatially distant places (or may 
be near and separated later) and we aim at extracting an 
entanglement between the two-cavity states by repeat- 
edly bringing the two-level atom into interaction with 
them and then selecting a particular state of the atom 
by measurements. 

The ideas of generating Q and of transferring Q en- 
tanglement in two-cavity system have already been pro- 
posed and studied in the context of the cavity quantum 
electrodynamics (CQED) [ljj, spectacularly developed 
over the last two decades both in the microwave ^l| and 
optical |l2j domains. Entanglement is generated from a 
properly prepared pure state in the former |8( and an 
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initially prepared entanglement in one cavity is trans- 
formed into another entanglement between the two cav- 
ities in the latter @, via successive interactions with a 
two-level atom. The atom plays the role of a "mediator" 
or "transformer" of entanglement. A similar, but more 
complicated role is sought for the atom in the present 
scheme, because its interactions with the cavities and 
the measurements of its state are expected to enable us 
to extract an entangled state, that is to produce an en- 
tanglement distillation, irrespectively of the initial states 
of the two cavities. 

For simplicity, suppose that the two cavities, A and B, 
are identical and their interaction with a two-level atom 
X is well described by the Jaynes-Cummings Hamilto- 
nian |l3j | . Let a and b indicate the annihilation operators 
of the modes of the two cavities A and B, respectively. 
The free and interaction Hamiltonians are 

H = — (1 + cr 3 ) + wa*a + u)tfb, (4.1a) 

h xa = 9A(o-+a + cr-a f ), H^ B = g B (o-+b + er_6 T ). 

(4.1b) 

A state where the atom X is in up (down) state and the 
modes a and b are in the nth and mth levels, respectively 
is denoted as |T(|), n, m) (n, m = 0, 1, 2, . . .). 

Since our purpose is to extract a pure state not in prod- 
uct form but entangled, it turns out that simple processes 
like (|3.2(l for the three-qubit system would not work. 
Indeed, because of the choice of the Jaynes-Cummings 
(rotating- wave) interactions l|4.1bj) . the number operator 
(1 + <73)/2 + a) a + wb commutes with any of H , H^ A , 
and i?xB ' an( i therefore any state of the two-cavity sys- 
tem of the form |n,0), which is a product state, is easily 
seen to be an eigenstate of the time evolution operator 
constructed analogously to 1)3. 4f> . (If the down state of X 
is measured, product states |0,m) are found to be eigen- 
states of the relevant time evolution operator.) Thus a 
process different from (|3.2|) would be necessary for our 
purpose. 

The above consideration would suggest that, with in- 
teraction given by Ij4.1b|) , it would be better to select after 
measurement a state of X different from the initial state. 
However, at the same time, we need a procedure that can 
be repeated many times within the present framework of 
the purification mechanism. Thus we choose a procedure 
that can be described schematically as 

preparation in ||) 

— > interaction between X and A for £a 



ID 



-HI) 




IT) 

FIG. 2: A two-level atom X, prepared in J.), is brought to 
interaction with cavity modes a and b in the two cavities A 
and B successively, and its state is measured after the inter- 
actions. Atoms that are found in state ||) will be sent back 
to the cavities in the reversed order. The state of atom X is 
again measured, and if it is found in ||), the whole process is 
repeated; other events are discarded. 

— > free evolution for ta 

— > interaction between X and B for t B 

— > free evolution for tb 

— > projection to |f) 

— > free evolution for tb 

— > interaction between X and B for t B 

— > free evolution for ta 

— > interaction between X and A for £a 

— > projection to (4-2) 

See also Fig. 

This is clearly a generalization of the purification pro- 
cess, "projection" — ► "time evolution" — > "projection." 
Indeed, in the above scheme, "time evolution" is not 
meant in the usual sense, that is described by a total 
Hamiltonian. It is instead interrupted by another pro- 
jection. However, the condition under which the purifi- 
cation mechanism does work is essentially the same as in 
the ordinary cases and all what we have to do here is to 
investigate the relevant evolution operator correspond- 
ing to the above process (14.2(1 . It would be important to 
notice that the above choice of the initial and projected 
states for system X is not arbitrary. In fact if it were pre- 
pared in the up state, the procedure analogous to that 
described before does not work. The vacuum state of the 
two cavities, which is a product state, would indeed turn 
out to be an eigenstate. 

The relevant evolution for the above process (|4.2(l is 
represented by products of the time-evolution and pro- 
jection operators and each of them is easily evaluated. 
The only non-trivial operators are 



oo 

e -i( Ho +H^ A )t A = e -i^u A e -(«+iMA ( cos ^n+i) | Tj n )(1,n\-i sin |T, n)(|, n + 1| 

- i sin ^i n+1) ||, n + 1)<T, n\ + e^- cos n) (I, n|) (4.3a) 
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and 



,m| — isin</5g +1 \'\,m}{l,m+ 1| 



e -i(ffo+^ B )*B = e ^at B ^2 e - l ( m+1 ^ tB (cos^ m+1) |T,m>(T 

- i sin (/4 m+1) II, m + 1) (T, m| + e^' B cos ||, m) (|, m| 



m=0 



where angles ip A and </? B m '' are defined as 

Pa ) = QA^Vn, ^4 m) = 5b£bV™- 
It is an elementary task to evaluate in this case the relevant evolution operator V c , analogously to (|2.5|l . 



\4 = (|| e _i ( ffo + ff xA) t A e -i-H"orA e -i(-H"o+-H"xB)*B e -i-H"oTB||^|| f 

and its explicit expression reads 



- i H o t b - i ( ffo + J?x b ) 4 b p - i-ffo U p - i ( ff o + H x A ) t / 



(4.3b) 



(4.4) 



II), (4.5) 



v; 



E « 

n,?n— 



-2i(n+m)ujT 



I ■ 2 («) 2 (m+1) , 2 («) • 2 (™)m \/ I 

(sin y^'cos ' + cos ^'sin ip^ ')\n,m)(n,m\ 

, . M ■ (in) (m) I . 1 1 > / i 

+ sin^ cos if A sin cosing \n + l,m — l){n,m\ 

, . (n) fn— 1) • (m+1) (m+1) i -. , i \ / i 

-I- sin <^ cos <£> A sm^g cos^g \n — 1, m + 1) (n, m\ 
I 



(4.6) 



where T = t a + ta + t B + t b . 



(n + m = 0) belongs to zero-eigenvalue of 14 and we need 



It is manifest from l|4.6|) that in the product Hilbert not consider it, for |0, 0) = 0. This is closely related to 
space of the two cavities there are sectors of V c within the choice of the initial state of X and is the reason why 
which the action of V c is closed. These invariant sectors we must prepare X in the down state ||). 



are characterized by the number n + m. We have (n + 
m+ 1) states {\n + m, 0), . . . , |0, n+m)} for the (n + m)th 
sector. Notice that the singlet state (vacuum state) |0, 0) 



Let us turn our attention first to the doublet subspace 
(ri + m = 1). The action of V c on this subspace is easily 
read from l|4.6|l as 




/ • 2 (1) 2 (1) 
-2iuT / Sm Pa C0S PB 



(1) • (1) (1) 

ysm ip A sm (p^ cos </? B 



(l) . (i) (i 
sm ip A sm </? B cos </? B 

• 2 (i) 
sm ip^' 




(4.7) 



Observe that the determinant of this matrix always 
vanishes, which means that one of the eigenvalues is 
zero and the other is given by the trace of the ma- 



trix 



-2iuT( 



(sin 2 (p" A cos 2 ip^' Pft')- Therefore we 

have a possibility of obtaining the largest (in magnitude) 
eigenvalue by adjusting the parameter ^a^a so that 



(i) 
sm ip A 



sm(g A tA) = ±1. 



(4.8) 



(The possibility cos<p B ^ = would result, not in an en- 
tanglement distillation, but in a product-state purifica- 
tion.) In such a case, the above eigenvalue equation is 




,-2iwT 




' (1) -r. ■ (1)\ 

cos <p B ' =F sm \ 

_i_ • (i) (i) j 

V ±smt/? B ' cost/? B / 

/ (i) _i_ • (i) 
/ cos ip B ± sm <p B 

\=F sin cos ip^ 



from which it is clear that the entangled state ) = 
cos ipjp |1, 0) ± sin |0, 1) can be extracted with the 
maximal probability by this setup. Notice that it still 
remains the freedom to adjust the value of y B = <7b£b- 
The remaining task is to check whether there are other 
eigenstates of V c belonging to eigenvalues with unit mag- 
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nitude, under the condition (|4.8(l . Consider the invariant 
sector characterized byfc = n + m>l that is composed 



V c 



where matrix elements Cj , dj can be read from l|4.6|) as 

Cj = sin 2 $ cos 2 <p<£- j+1) + cos 2 p^ sin 2 , 

(4.11) 

dj = sin^ costpV-V sin^' +1) cos^' +1) . (4.12) 

It is important to notice that the condition (|4.8|l im- 
plies that the element d 2 = sin tp^ cos ipj^ sin ip^ ^ x 

cos ipg* ^ vanishes, irrespectively of p^ ^ and thus the 
sector further splits into two subsectors {\k, 0), . . . , |2, A:— 
2)} and {|1, k — 1), |0, fc)}. Furthermore, it is easily seen 
that the entangled state in the latter subspace of the form 
l^fc ) = cos ip^ |1, k — 1) ± sinip^\0,k) has the eigen- 
value _ e - 2lfew " ; while, as shown in the Appendix, no 
eigenstate in the former subspace {\k, 0), . . . , |2, k — 2)} 
belongs to a unit (in magnitude) eigenvalue (if k is 
smaller than 9). 

We have seen that there are, for any fc-sector, many 
entangled states ) = cos^g |1, fc— 1) isin^g |0, k) 
{k = 1,2,...) (that increase, in general, in higher k- 
sectors) extracted with the optimal probabilities by the 
process 1)4. 2|) . Repeated interactions of the two-level 
atom X in the cavities A and B and the prescribed mea- 
surements (projections) certainly bring us with a statis- 
tical (classical) mixture of these entangled states. The 
situation would not be considered completely satisfac- 
tory, since we would not be able to distill a single en- 
tangled state by the process l|4.2|) . There is, however, 
a way out of this difficulty. We may prepare such an 
initial (mixed) state of A+B that contains only those 
sectors with relatively small fcs. Such a preparation of 
the initial state would effectively eliminate the possibil- 
ity of obtaining other states than \$>^) after performing 
the process (|4.2f) . For example, we may consider the 
following preparation procedure, which is nothing but 
a purification process applied to cavity B. We send a 
two-level atom prepared in down state ||) to cavity B. 
After its interaction, which is again assumed to be of the 
form (|4.1b|) . with B, the atom is measured and only those 
events in which it is found in the state ||) are retained. 
This process is to be repeated many times and the result- 
ing state of A+B, which will be used as the initial state 



of k+ 1 states {|fc, 0), . . . , |0, k)}. The action of V c on this 
sector is represented by the following matrix [see l|4.6l) ] 



(4.10) 



I 

for the following entanglement distillation process (|4.2|l . 
would be dominated by the state pab (0) ~ Pa (0) ® |0) (0| , 
since the vacuum state of system B is the unique eigen- 
state of the relevant evolution operator ~ (||e _lffxB * |J.) 
belonging to eigenvalue unity if no fine tunings are made 
on the parameters. After having prepared the state 
Pab(0), we repeat the process (|4.2|) under the condition 
(|4.8|) . We would finally end up with the single entan- 
gled state l^c 1 ^), because our initial state pab(0) satisfies 
(^ ) |pA B (0)|^ fe) ) = (^ fe) |pAB(Q)|^ fe) ) = for k > 1, 
where (<E>c^ | is the left-eigenstate corresponding to \^c) • 
Concluding this section we wish to give typical values 
of relevant parameters under the aspect of the possibility 
of implementing our proposal in laboratory. We concen- 
trate on the estimation of the total duration T of the 
experiment. To this end, we choose to be in the context 
of microwave CQED where both the geometrical arrange- 
ment of the experimental set up and the intensity of the 
atom-field coupling regime, seem more favorable to our 
proposal. Let us first note that the typical atom-field 
coupling constant g (<?a or or) can be chosen in such a 
way that g ~ 10 4 -10 5 s -1 [Tj|- Moreover the lifetime of 
a Rydberg atom is > 10 _2 s ^5[. As for the quality 
factor Q of the cavities currently used in labora tory , we 
quote typical values of the order of 10 8 -10 10 |14lll5l| cor- 
responding to a cavity damping time 1 ms-1 s. Consider- 
ing that in our case £a ~ £b ~ g l , the total duration T 
of the experiment may be estimated and turns out to be 
compatible with the entanglement distillation proposed 
in this section. 



V. SUMMARY 

In this paper, the idea of extracting entangled states 
among systems located at spatially separated places, ir- 
respectively of their initial states, has been proposed and 
applied to simple systems to show the potentiality of a 
novel measurement- based purification scheme .3]. The 
establishment of entanglement distillation relies on the 
successive interactions between the systems under con- 
sideration and the so-called "mediator" quantum system. 



/ IM) \ 
lfe-i,i' 



\2,k-2) 
|l,fc-l) 
|0,fc) 



-2ikuT 



/ Cfe d k 
dk Ck-i 






Vo 



0' 




c 2 d 2 
d 2 ci di 
di c / 



/ IM) \ 
I* -1,1 



\2,k-2) 
\l,k-l) 

V |o,fc) J 
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In the first example, it is demonstrated that the entangle- 
ment between the two qubit states is possible via their in- 
teractions with another qubit, which plays the role of the 
entanglement mediator, with the maximal yield (optimal 
distillation) . In the case of distillation of cavity- mode en- 
tanglement, however, a modification of the original sim- 
ple scheme, that is, "interaction — > measurement — > in- 
teraction is required and the modified procedure 
(|4.2|) turns out to result in the entangled state |vE'c 1 ''), 
after an appropriate preparation of the initial state. We 
stress that there would also be a possibility of obtaining 
an entangled state not only in the lowest sector k = 1 
but also in the higher k sector, provided an appropriate 
initial state be prepared. 

It would be worth stressing that in spite of such a mod- 
ification required in the second example, the underlying 
notion is still the same: the action of a measurement (rep- 
resented by a projection operator, for simplicity) causes 
an essential and critical dynamical change, not only in 
the system measured, but also in the others interacting 
with the former. Since the notion is so general, one can 
devise various applications of this scheme in many differ- 
ent situations. The examples explored in this paper are 
just two of them and further applications will be reported 
elsewhere. 

Finally we add some comments on the practical setup 
of our proposal. In both schemes reported in the previous 
sections, the entanglement mediator is an atom appro- 
priately prepared before entering in interaction with the 
two subsystems to be entangled and subjected to a condi- 
tional measurement of its internal state at the end of the 
two successive coupling processes. We wish to stress that 
to realize in practice the required "many crosses" scheme 
one may synchronize the injection of the jth atom of the 
sequence into the process with the successful outcome 
of the internal state measurement of the (j — l)th one. 
This kind of experimental setup might be preferable to 
the conceptually simpler one based on the idea of using 
always the same atom reversing its direction of motion 
at exit to re-inject it into the process. Thus the repre- 
sentations in Figs. 1 and 2, as far as the aspect under 
scrutiny is concerned, have been reported only for the 
sake of simplicity. 
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APPENDIX 



2) 



In this Appendix, a symmetric matrix of the form [k > 



/ Cfc d k 
dk Cfc-i dk-i 
d k -\ c fc _2 





Vo 



... 




... 





... 





... c 3 


d 3 


...da 


J 



(A.l) 



with the matrix elements [see l|4.11fl and (|4.12|) ] 



• 2 U) 2 (k-j+1) 

sin tpy cos ip B 



cos 2 ^sin 2 <pg- j \ 

(A.2a) 



7 • {j) • (k-j+l) {k-j+l) 

dj = sin ip£ cos (pY ' sin (p B cos <p B , 

(A.2b) 

is investigated with particular attention on its eigenstates 
belonging to the eigenvalues with unit magnitudes. Since 
this matrix is real and symmetric, its eigenvalues are all 
real, and the eigenvalues of relevance in the framework 
of our procedure here are ±1. 

For the first possibility +1, let us consider the deter- 
minant U (fc > i > 2) defined by 



(H - 1 di 
di Cj_i - 1 di^i 
di^i Ci_ 2 - 1 



C3 - 1 d 3 

d 3 c 2 - 1 



(A.3) 

It is easy to see that the particular form of Ii and the 
definitions of Cj and dj in (| A. 2|> lead to a recursion rela- 
tion 

7-1 2 (i) 2 (fc— i) t 
I i + COS ipy COS (p B 'Ii-i 

• 2 (i) • 2 (fc-i+1) 

= — sin ip A sin ip B 



I r 2 (i— 1) 2 (k—i+l) t 

X I + COS Z (f\ ' COS" 2 ^ ' Ji_ 2 



for k > i > 4. This is further reduced to 



(A.4) 



T 2 (i) 2 (fc— i) t 

Ii = -cos ipy cos ip B 'Ii-i 



i-l 



,2 ,„0+i) -2 ,„(fc-i) 



(-If 1 JJsinVx ' sm"^ 

(fc>»>2). (A.5) 
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If we set It = (— l) l+1 Pi, then Pi is found to be positive This relation is easily solved to yield the explicit form of 
scmi-dcfinitc and to satisfy If. = (— l) fc+1 P/ c with 

P = cos 2 $ cos 2 yg-^Pi^+Jl sin 2 sin 2 <pg-» . 

3 = 1 

(A.6) 

I 



n 2 (2) 2 (fc) 2 (1) 2 (fc-2) 

P fe = COS ip A ■ ■ ■ COS <f\ ' COS V?B ' COS cpg ' 

fc-1 

, • 2 (2) -2 (n) 2 (n+1) 2 (fe) 2 (1) 

+ > sin sm cos cos Va cos 

n=2 

, • 2 (2) • 2 (fc) • 2 (1) -2 (fc-1) 

+ sm ip A ■ ■ ■ sm y> A sm ip^ ■ ■ ■ sm y> B 



2 (fc-n-l) • 2 (fc-™+l) . 2 (fc-i) 
■ cos sm • • • sm <y9 B 



(A.7) 



Since each term in (|A.7|I has the same sign to each 
other, vanishing of Ik , which is nothing but the condition 
for the matrix (|A.1|I to possess eigenstates that belong to 
eigenvalue unity, is equivalent to that of each term. This 
means that there are k conditions for three parameters 
Sa^a, 3b^b and k and it seems impossible to have a van- 
ishing Ik in general, unless most of the conditions are 
simultaneously satisfied. If we choose a particular value 
for gA^A, say g\tA = tt/2 as in l|4.8|) . however, all Ik with 
k > 9 vanish because each term in (|A.7|) contains either 

sin 2 ip)^ or cos 2 ip A \ each of which vanishes. 

On the other hand, the matrix (|A.1|I is shown to have 
no eigenstate belonging to the eigenvalue —1. Consider 
the following determinant 



> CiJi-l — d; Jj_2 



a + 1 
di 




di 
Cj_i + 1 di-i 

di-l Ci-2 + 1 



where Cj and dj 







c 3 + l 
da 







d 3 
c 2 + 1 



are again given in (| A. 2|> . Since 



(A.8) 



> 



0, d'j < 1 and therefore J 2 — c-i + 1 > and J3 = (03 + 
l)(c 2 + 1) — d\ > 0, let us assume that all J/s are positive 



definite for I up to i — 

Ji = (Cj + 1)J<-1 



1. Then it follows that 
— d]ji^2 



■ 2 (i) 2 (fc-'i+l) 

= sm </? A cos yjg 

/ 7 2 (i-l) • 2 (fc—i+1) 7 

x ( Ji-i - cos ip\ sm tpg 'Ji_ 2 

2 (») • 2 (fc-») 7 

+ cos ^ sin y>g Ji-x- 
This relation recursively yields inequalities 

T 2 (i) • 2 (fc— i) 7 

J, - cos <p A sm (^g ; Jj_i 

. • 2 K) 2 (fc-i'+l) 
> Sin lf A ' COS Iy5g 

/ 7 2 («-l) • 2 (fc-i+1) 7 

x ( Jj_i - cos z ip A ' sm z y> B ' Ji_ 2 



(A.9) 



. . 2 (*) 2 (fc-i+i) 
> sm ip A cos (^g ' 



• 2 (4) 2 (fc-3) 

■ sm ip A ' cos ^g ' 



( 7 2 (3) ■ 2 (fc-3) 7 

x I J3 — cos </3 A sm <^ B J 2 



(A.10) 



The last factor on the right most hand is shown to be 
positive definite 



7 2 (3) . 2 (fc-3) 7 

J 3 - cos ip A sm 93 B J 2 

> (sin 2 tp { l } cos 2 (p B fe ~ 2) + 1 ) J 2 - d'i > 0, (A. 11) 



which means that the quantity on the left-hand side in 
(|A.10|) is positive definite. We conclude that Ji is pos- 
itive definite and therefore Jk does not vanish, which 
completes the proof. 
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